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ABSTRACT 



The bright gamma-ray burst GRB050525a has been detected with the Swift 
observatory, providing unique multiwavelength coverage from the very earliest 
phases of the burst. The X-ray and optical/UV afterglow decay light curves both 
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exhibit a steeper slope ~ 0.15 days after the burst, indicative of a jet break. This 
jet break time combined with the total gamma-ray energy of the burst constrains 
the opening angle of the jet to be 3.2°. We derive an empirical 'time-lag' redshift 
from the BAT data of z — 0.69 ± 0.02, in good agreement with the spectroscopic 
redshift of 0.61. 

Prior to the jet break, the X-ray data can be modelled by a simple power law 
with index a = —1.2. However after 300 s the X-ray flux brightens by about 
30% compared to the power-law fit. The optical/UV data have a more complex 
decay, with evidence of a rapidly falling reverse shock component that dominates 
in the first minute or so, giving way to a flatter forward shock component at later 
times. The mult i wavelength X-ray/UV/Optical spectrum of the afterglow shows 
evidence for migration of the electron cooling frequency through the optical range 
within 25000s. The measured temporal decay and spectral indices in the X-ray 
and optical/UV regimes compare favourably with the standard fireball model 
for Gamma-ray bursts assuming expansion into a constant density interstellar 
medium. 

Subject headings: astrometry - galaxies: distances and redshifts - gamma rays: 
bursts - shock waves - X-rays: individual (GRB 050525a) 

1. Introduction 

The observation of Gamma-ray bursts has entered a new era with the launch of the 
Swift observatory (Gehrels et al. 2004). Swift provides rapid notification of GRB triggers to 
the ground using its sensitive Burst Alert Telescope (BAT; Barthelmy et al. 2005), and can 
make panchromatic observations of the burst and its afterglow by bringing its narrow field 
X-ray Telescope (XRT; Burrows et al. 2005a) and UltraViolet/Optical Telescope (UVOT; 
Roming et al. 2005a) to bear within about 1 minute of the burst going off. 

In this paper we describe Swift observations of GRB050525a, the first bright, low- 
redshift burst to have been observed with Swift with both its narrow field instruments. This 
has resulted in the most complete multiwavelength description of the early evolution of a 
burst afterglow yet obtained. We describe the light curve decay in the X-ray and optical/UV 
bands, and compare these with the predictions of theoretical models. We also consider how 
the multiwavelength spectrum of the burst afterglow evolves with time, making use of the 
simultaneous X-ray/UV/optical coverage. We combine the measurements of the afterglow 
properties with the well constrained measurements of the prompt Gamma-ray emission to 
constrain the geometry of the burst emission. We present an analysis of the prompt gamma- 
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ray emission in section 2.1. This is followed by an analysis of the time evolution and spectrum 
of the X-ray data in section 2.2, adding to this the optical and ultraviolet data in section 2.3 
to build up the multiwavelength picture. The results are discussed in section 3. 

2. Observations and Analysis 
2.1. BAT 

At 00:02:53.26 UT, the Swift Burst Alert Telescope (BAT) triggered and located on- 
board GRB050525a (trigger #130088). The BAT location calculated on-board was RA, DEC 
278.144, +26.340 (J2000) with an uncertainty of 4 arcmin (radius, 3a, including estimated 
systematic uncertainty). This was later refined (see Figure 1) to 278.140, +26.344 ± 0.5' 
(95% containment, the position corresponding to 18h 32m 34s, +26d 20' 38") using the full 
data set downloaded through the Malindi groundstation ~1 hour after the burst. The refined 
position is 26 arcsec from the on-board location and 21 arcsec from the UVOT location. The 
burst was 28° off the BAT bore sight (85% coded) and was 148cx in the image-domain full data 
set. The spacecraft completed the automated slew to the burst at 00:04:58 UT (75 sec after 
the trigger), and XRT and UVOT began their standard follow-up observation sequences. 

Figure 2 shows the BAT light curves in the 4 standard BAT energy bands. The burst 
was very bright and was detected in all 4 energy bands with a peak rate of ~80,000 cts s -1 . 
The burst has two main peaks with several smaller peaks within the main structure. The lag 
between the lightcurves in the 100-350 keV and the 25-50 keV bands is r = 0.124 ± 0.006 s 
(the high energy band leads the low energy band). The total duration was ~12 s (T90 is 
8.8+0.5 s). The successful trigger criterion for this burst was in the 25-100 keV band with 
a duration of 64 msec. 

Another very bright burst occurred ~2800 s later (GRB050525b; Golenetskii et al. 
2005). This was outside the field-of-view (FOV) of INTEGRAL'S IBIS instrument whereas 
GRB050525a was inside the IBIS FOV. Thus we can be confident that this second burst was 
not associated with GRB050525a (Mereghetti et al. 2005). 

We fit BAT spectra using the current analysis techniques developed by the instrument 
team. Systematic effects are still being resolved this early in the mission, and consequently 
the fits described here should be regarded as preliminary. To reproduce the known spectrum 
of the Crab, empirical correction factors have been applied to the pre-launch response matrix; 
this analysis also resulted in an estimate of the energy-dependent systematic error that should 
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be added to the statistical error. The corrections and the systematic error are posted. 1 The 
count spectra are binned between 16 and 148.8 keV in ~2 keV bins. Table 1 summarizes the 
spectral fits to the entire burst (12.80 s) and to the peak 1 s, with 90% confidence limits. We 
fit the count spectra with three nested models, here presented as energy spectra 2 : a power 
law F(E) oc E 13 ; a power law with an exponential cutoff F(E) oc E@ exp[— E/Eq\; and the 
'Band' model (Band et al. 1993), a low energy power law with an exponential cutoff that 
transitions into a high energy power law F(E) oc E^ 2 . The peak energy E p = (1 + (3)E is 
both physically more relevant and less correlated with (3 than E ; E p is the energy of the 
peak of E F(E) oc vf v . The power law with an exponential cutoff is the same as the Band 
model with (3 2 = — oo, and the power law model is the same as the other two models with 
E = oo. 

As can be seen from these fits, the power law fit is definitely inferior to the other two 
spectral models, and the second and third spectral models are equivalent in this case since 
the transition to the high energy power law at Ef, = {(5 — ^Eq = {(5 — (32)E P / ((3 + 1) is above 
150 keV (the energy of the highest PHA channel that is fitted); indeed, for the observed (3 
and E p values, the fits are insensitive to (3 2 < —2.5. 

Konus-Wind (Golenetskii et al. 2005) and INTEGRAL (Gotz et al. 2005) both detected 
GRB050525a. Golenetskii et al. (2005) fit the Konus-Wind 20 keV to 1 MeV spectrum for 
the entire burst with a power law with an exponential cutoff model and find (3 = — 0.10±0.05, 
and E p = 84.1 ± 1.7 keV, consistent with the Swift fit. 

The spectral fits can be integrated over energy and time to give a fluence S( 15-350 keV) 
= (2.01 ± 0.05) x 10~ 5 erg cm~ 2 and a peak flux of P = (15-350 keV) = 47.7 ± 1.2 photons 
cm~ 2 s _1 . Using the observed redshift of z — 0.606 (Foley et al. 2005), the spectral fits and 
a standard cosmology (f_ m =0.3, Q^=0.7, and H = 70 km s -1 Mpc -1 ), we find E iso = 2.3 x 
10 52 erg and L iso = 7.8 x 10 51 erg s" 1 . These two derived quantities are the total energy and 
peak luminosity that would have been radiated if the observed flux were emitted isotropically. 
As will be discussed below (section 3) in considering the possible jet break, the emission was 
most likely beamed and the actual emitted energy and luminosity were consequently smaller. 
Using the redshift, the peak energies in the burst frame are E pp = 164.5 ± 11.6 keV for the 
peak of the lightcurve and E pt = 126.6 ± 5.5 keV for the entire burst. 



1 http://legacy. gsfc.nasa.gov/docs/swift/analysis/bat_digest. html 

2 We use the notation (3 here and throughout the paper to denote the power law index of the energy 
spectrum of the source. 
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Table 1. BAT Spectral Fits 





Entire Burst 


Peak Flux 




Parameter 


Power Law a Power Law, Band Model 


Power Law a Power Law, Band Model c 






Cutoff b 


Cutoff b 
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n m +011 
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-0.42 c 


n 45+ - 14 


-0 45+ 14 
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-7.84l?;S 






o 97+7.14 
°- z ' -0.73 


E p z 
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102.4+^ 
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i5.otl:2i h 


19.35° fe 


7 62+ 3 91 1 


44.00±|;§g h 


X 2 /dof 


181.7/57 


15.2/56 


15.2/55 


169.8/57 


31.6/56 


31.6/55 



a Power law model, F(E) oc E 13 . 

b Power law with an exponential cutoff, F(E) oc E" exp[— E/Eq\. 

c Band model (Band et al. 1993), a low energy power law with an exponential cutoff transitioning to a high 
energy power law F(E) oc E^ . 

d The low energy spectral index. 

e Fit too poor to produce uncertainty range. 

f The high energy spectral index. The fit is insensitive to fii < —2.5 for the fitted E p . 

«The energy of the peak of EN(E) oc vf v , and E p = (2 + /3)E . 

h The normalization of the spectrum at 50 keV, in keV cm~ 2 s _1 keV -1 . 

'The normalization of the spectrum at 1 keV, in keV cm -2 s _1 keV -1 . 
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2.2. XRT data 

The XRT observations started on 2005 May 25 at 00:04:08 UT, 75 seconds after the trig- 
ger. The XRT was undergoing engineering tests at the time of the burst and the Auto State 
(which selects the observing mode based on source brightness) was disabled. A sequence of 
3 frames in Photon Counting (PC) mode were taken first, between T+75s to T+83s. These 
early PC mode data suffer from severe pile-up because of the brightness of the source, and 
have been excluded from the analysis. 

Shortly afterward the instrument was put into Auto State and a 2.5 s exposure in Image 
mode was taken starting at 00:04:58 UT, 125 seconds after the BAT burst trigger. A bright 
X-ray source was found near the center of the field of view. The refined X-ray coordinates 
are RA(J2000) = 18 h 32 m 32 s .6 and Dec(J2000) = 26°20'18", with an estimated positional 
uncertainty of 6 arc seconds (90% confidence level). The XRT coordinates are 4 arc seconds 
from the UVOT position of the optical counterpart (see below). 

Following the automatic sequence of readout modes, designed to avoid any pile-up effect, 
the instrument was configured in the Photodiode (PD) readout mode starting from 00:05:01 
UT (T+128 s). Due to the engineering tests mentioned above, the instrument remained in 
PD mode until 00:20:21 UT (T+1048 s) even when the afterglow brightness was below the 
nominal flux threshold used for this mode. 

Data in PC mode, the most sensitive XRT operational mode, were taken starting at 
01:40:32 UT (T+5859 s) until 11:39:11 UT (T+41778 s). The GRB050525a field was re- 
observed in PC mode on several occasions up to late June 2005. 

2.2.1. Temporal analysis 

For the PD mode, events for the temporal analysis were selected in the 0.4-4.5 keV en- 
ergy band to avoid contamination from the calibration sources. The standard grade selection 
for this mode (0-5) was used. The count rate was converted to unabsorbed 2-10 keV flux 
using the best fit spectral model (see below). 

For PC data, events were selected in the 0.3-10 keV band and grades 0-12 were used in 
the analysis. Furthermore, since 2-D spatial information is available for this mode, photons 
within a circle of 10 pixel (~24 arcsec) radius, which encloses about 80% of the PSF at 1.5 
keV, were extracted, centered on the source position. The background was estimated from a 
nearby source-free circular region with 50 pixel radius. Again, the count rate in the 0.3-10 
keV band was converted to unabsorbed 2-10 keV flux using the best fit spectral model. 
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a time since trigger 

b flux in units of 10~ n erg cm~ 2 s _1 2-10 keV 
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Table 3. UVOT multicolor data 
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The background subtracted 2-10 keV light curve in the time interval T+128 s - T+1048 s 
(PD mode) is shown in Figure 3 (inset). The X-ray afterglow of GRB 050525 is clearly fading. 
The early afterglow decay was first fitted with a single power-law model, resulting in a best 
fit decay index a = —0.95 ± 0.03, with xl — 1-17 (42 dof). Inspection of the residuals to 
the best fit model suggests that a flattening of the decay curve or a re-brightening of the 
source occurs at ~300 seconds after the trigger. A better fit is provided by a broken power 
law model with slopes a 1: a 2 and a break at t b . This model gave xl — 0-98 (40 dof), with 
best fit parameters «i = —1.231^02) a z = — 0-91 and t b = 203 s. 

Again, however, the residuals suggest sytematic deviations from this model. We thus 
tried a broken power law with two temporal breaks. This model provided a very good fit to 
the data, with xl — 0-72 (38 dof) and is plotted in Figure 3 (inset) as a solid line. The best 
fit parameters are cti = —1.19, t\ = 282 s, a 2 = —0.30, t\ = 359 s, and a 3 = —1.02. 

Next, we fitted the X-ray data taken in PC mode at times more than 5000 s after 
the trigger. We first used a single power-law model, obtaining a best fit decay index a = 
— 1.51 ± 0.07, with xl — 1-40 (12 dof). The poor fit is the result of a clear steepening of 
the light curve with time. We thus tried a broken power law model. The model provided a 
very good fit with xl — 0-97 (10 dof) and best fit parameters a± = —1.16, a 2 = —1.62 and 
t b = 13177 s. 

Finally, we tried fitting the total light curve derived from the combined PD and PC 
mode data (see Figure 3). We find that the power law fit to the pre-brightening PD mode 
data (T<280 s) extrapolates well to the pre-break PC mode data. Moreover the decay index 
before 280 s agrees well with that of the PC mode data before the 13ks break. In contrast, 
if we extrapolate the post brightening PD mode data to later times using the best fit slope, 
a significant excess is predicted compared with the measured PC mode data. To join the 
post brightening PD mode data to the PC mode data requires a model with at least two 
temporal breaks, which are not constrained because of the intervening gap in X-ray coverage. 
We conclude that the brightening at about 280 s in the PD mode data represents a flare in 
the X-ray flux, possibly similar to the sometimes much larger flares that are seen at early 
times in other bursts (Burrows et al 2005b; Piro et al. 2005), and that the flux returns to 
the pre-flare decay curve prior to the start of our PC mode data. 

We thus fit the combined PD and PC mode data excluding PD data at times t >T+288 s 
(green points in Figure 3). A broken power law model provided a good fit (solid line of 
Figure 3), with with xl — 0.50 (25 dof) and best fit parameters a\ = —1.20 ± 0.03, a 2 = 
-1.62^6 and t b = 137261^ s. The break time is thus ~ 3.8 hours. 

The complete XRT data are recorded in Table 2. 
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Table 3 — Continued 



T(mid) a T(cxp) 
(s) (s) 


Mag 




Flux Density* 3 


T(mid) a 


V v 1~V 


Mag 


Flux Density 13 


468 


10. 


14.71 


± 


0. 


.19 


43.5 ± 8.3 


28664 


194. 


18.85 ± 0.30 


1.0 ± 0.3 


552 


10. 


14.97 




0. 


.20 


34.3 ± 6.9 


29380 


900. 


18.79 ± 0.20 


1.0 ± 0.2 


637 


10. 


15.01 


± 


0. 


.20 


33.0 ± 6.7 


41040 


758. 


19.99 ± 0.54 


0.3 ± 0.2 


721 


10. 


15.14 


± 


0. 


.20 


29.3 ± 5.9 


52994 


780. 


> 19.68 


< 0.4 


805 


10. 


15.46 


± 


0. 


.21 


21.8 ± 4.7 


128928 


1548. 


> 20.05 


< 0.3 


UVW1 (260 


nm) filter 


















202 


10. 


13.61 


± 


0. 


.07 


162.8 ± 10.4 


877 


10. 


16.15 ± 0.24 


15.7 ± 3.5 


286 


10. 


14.18 


± 


0. 


09 


96.3 ± 8.0 


961 


10. 


15.76 ± 0.20 


22.6 ± 4.1 


370 


10. 


14.77 


± 


0. 


.12 


56.2 ± 6.2 


6761 


638. 


17.33 ± 0.06 


5.3 ± 0.3 


455 


10. 


15.10 


± 


0. 


.14 


41.5 ± 5.3 


28111 


900. 


19.49 ± 0.19 


0.7 ± 0.1 


539 


10. 


15.22 


-1- 


o 

\J ■ 


.15 


37.0 ± 5.1 


40585 


900. 


20.33 ± 0.39 


0.3 ± 0.1 


624 


10. 


15.26 


± 


0. 


.15 


35.6 ± 4.9 


52100 


898. 


> 21.05 


< 0.3 


708 


10. 


15.22 


± 


0. 


.15 


36.9 ± 5.1 


63738 


900. 


> 21.20 


< 0.3 


792 


10. 


15.88 


± 


0. 


21 


20.1 ± 3.8 










UVM2 (220 


nm) filter 


















187 


10. 


13.54 


± 


0. 


10 


218.9 ± 20.0 


863 


10. 


15.57 ± 0.26 


33.8 ± 8.1 


272 


10. 


14.32 


± 


0. 


.14 


106.8 ± 14.1 


947 


10. 


15.66 ± 0.27 


31.2 ± 7.6 


356 


10. 


14.60 


± 


0. 


.16 


82.4 ± 12.4 


5984 


900. 


17.10 ± 0.06 


8.2 ± 0.4 


440 


10. 


14.94 


± 


0. 


19 


60.3 ± 10.7 


24007 


864. 


18.96 ± 0.19 


1.5 ± 0.3 


525 


10. 


15.84 


± 


0. 


.30 


26.4 ± 7.2 


39683 


897. 


19.72 ± 0.30 


0.7 ± 0.2 


609 


10. 


14.90 


± 


0. 


19 


62.4 ± 10.8 


51257 


900. 


20.34 ± 0.48 


0.4 ± 0.2 


694 


10. 


15.18 


± 


0. 


21 


48.3 ± 9.4 


74798 


931. 


> 21.59 


< 0.6 


778 


10. 


15.51 


± 


0. 


.26 


35.9 ± 8.5 










UVW2 (198 


nm) filter 


















245 


10. 


14.69 


± 


0. 


.11 


88.4 ± 9.5 


835 


10. 


16.43 ± 0.27 


17.8 ± 4.5 


329 


10. 


15.30 


± 


0. 


.16 


50.2 ± 7.2 


920 


10. 


16.21 ± 0.25 


21.8 ± 4.9 


413 


10. 


15.13 


± 


0. 


14 


58.6 ± 7.7 


1140 


100. 


16.35 ± 0.08 


19.1 ± 1.5 


498 


10. 


15.60 


± 


0. 


18 


38.2 ± 6.4 


17973 


882. 


18.99 ± 0.13 


1.7 ± 0.2 


582 


10. 


15.22 


± 


0. 


.15 


54.3 ± 7.5 


34806 


900. 


20.05 ± 0.25 


0.6 ± 0.1 


666 


10. 


15.80 


± 


0. 


20 


31.7 ± 5.8 


52162 


1800. 


> 21.76 


< 0.3 


751 


10. 


15.72 


± 


0. 


19 


34.3 ± 5.9 











a time since trigger 

b flux in units of 10 -16 erg cm -2 s _1 A -1 
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2.2.2. Spectral analysis 

The Photodiode spectra in the time intervals T+128 s - T+288 s ( "pre- flare" ) and 
T+288s - T+1048s ("flare") have been extracted in the 0.4-4.5 keV band. The spectra were 
binned to ensure a minimum of 20 counts per bin and were fitted using the XSPEC package 
(v. 11.3.2; Arnaud et al. 1996). 

An absorbed power law model fits both spectra well. For the "pre-flare" spectrum the 
best fit hydrogen-equivalent column density is Ah = 2.00 ± 0.25 x 10 21 cm~ 2 and the energy 
index is (3 = -0.99 ± 0.09 (x r 2 = 0.94, 176 dof). The best fit parameters for the "flare" 
spectrum are A H = 2.07 ± 0.18 x 10 21 cm" 2 , (3 = -0.96 ± 0.07 = 1.04, 238 dof). 

Freezing Ah to the Galactic value (9.0 x 10 20 cm -2 ) and adding a hydrogen-equivalent 
column density at the redshift of the GRB (z = 0.606) we found A H = 2.2 ±0.5 x 10 21 cm -2 , 
P = -0.90 ± 0.08 (x r 2 = 0.96, 176 dof) for the "pre-flare" spectrum and A H = 2.5 ± 0.4 x 
10 21 cm" 2 , = -0.89 ± 0.06 (x r 2 = 1.01, 238 dof) for the "flare" spectrum. 

Events in PC mode from the time interval T+5859 s to T+41778 s were extracted from 
the same circular region used in the temporal analysis. A further selection on XRT event 
grades 0-4 (i.e. single and double pixel events) was applied to the data. Again, the spectrum 
was binned to ensure a minimum of 20 counts per bin. Energy channels below 0.3 keV and 
above 10.0 keV were excluded. 

An absorbed power law model fits the PC data well (Xr — 0.70, 21 dof), with A H = 
1.5lg;? x 10 21 cm" 2 and p = -1.10l°;H 

No evidence is found for evolution of the X-ray spectral shape with time (cf. Chincarini 
et al. 2005). 



2.3. UVOT data 

The Swift UVOT began taking data in the V filter ~ 65 seconds after the burst trigger, 
immediately after the GRB came into its field of view. After about ten seconds the spacecraft 
attitude had settled and UVOT began a 100s 'finding chart' exposure in V. Thereafter the 
instrument cycled through each of six color filters, V, B and U together with filters defining 
three ultraviolet passbands, UVW1, UVM2 & UVW2 with central wavelengths of 260nm, 
220nm & 193nm respectively. The exposure duration per filter was initially 10 seconds, 
subsequently increasing to 100s and then 900s at predetermined times after the trigger. Data 
were taken in both 'event mode', in which the time and detector position of each individual 
photon is recorded, and in 'image mode', where the image is accumulated onboard, discarding 
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Fig. 1. — This Figure shows the field containing GRB 050525 A, marking the location of 
the optical afterglow (the object picked out by the arrow). The BAT error circle is based 
on the refined ground-based analysis (Cummings et al. 2005), and this has a radius of 0.5 
arcminutes. The XRT error circle has a radius of 6 arcseconds and the central position has 
been revised compared to that reported in Band et al. (2005) after ground analysis. 
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Fig. 2.— BAT lightcurve of GRB050525a in the 15-25 keV (solid), 25-50 keV (dots), 50-100 
keV (short dashes), and 100-350 keV (long dashes) bands. 
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Fig. 3. — XRT decay light curve of GRB050525a including both Photodiode Mode 
(T<2000 s) and Photon Counting Mode (T>2000 s) data. The solid line is a broken power 
law fit to the combined data excluding those Photodiode Mode points colored green (see 
text). The dash-dot line is shown for illustration and has a slope of a = —2.2, which is the 
value expected from simple modelling of a jet break (see text). The lower inset shows the 
data taken in photodiode mode only, during the first ~1000 s after the BAT trigger. The 
solid line is a fit to the data with a power law model that includes two temporal breaks to 
different decay rates. The dashed line is an extrapolation of a simple power law fit (single 
slope) to the first segment of data prior to about 300 s. The upper inset shows the residuals 
with respect to the two power-law model fit to all the data, expressed as a percentage of the 
predicted model flux. 
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Seconds since 2005 May 25.0020023 UT 



Fig. 4. — This Figure shows the photometry (solid circles) in the six UVOT broadband 
filters. Arrows indicate 3cr upper limits to the photometry. A model light curve (red line) 
with oil = -1.56, a 2 = -0.62, ay = -1.76, t t = 432 s (7.2 min), and t d = 13133 s (0.152 
days) is shown. The green line indicates a rebrightening model with a\ = —2.14, a 2 = —1.04, 
a>j = —1.73, t t = 39 s, and tj = 17 747 s (0.205 days). The flux has been normalized for each 
filter. 



5 




Fig. 5. — Comparison of flux decay in the X-ray and UVOT bands. The UVOT data have 
been normalised in the interval up to T+ 1000s, and the data taken through different filters 
are distinguished by color. The relative normalisation of the X-ray and optical/UV data is 
arbitrary. The best fit broken power law model is plotted through the X-ray data. The best 
fit double power law with break is plotted through the UVOT data (see text). The dashed 
line has the same post break slope as the X-ray data. The dotted line is the best fit model 
with a constant flux added corresponding to the value measured by Soderberg et al. (2005) 
using HST/ACS (see text). 
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Seconds since 2005 May 25,0020023 UT 



Fig. 6. — This Figure shows the UVOT photometry (solid circles) normalized so that the 
early time flux is the same in each filter. The red and green lines are the models from 
Figure 4. The yellow line is the best-fitting single power law. The open circles are the Klotz 
et al. (2005) data scaled to the same flux scale. The solid square is the HST/ ACS F625W 
data of Soderberg et al. (2005). Note that the HST / ACS data point is approximately three 
mag brighter than the predicted afterglow magnitude. 
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Frequency (Hz) 

Fig. 7.— UVOT and XRT data interpolated to the epochs T+250s, T+800s and T+25,000s, 
and best spectral fit models (solid lines). The dotted lines represent the intrinsic continuum 
of the source, before extinction and absorption from gas and dust in both the Milky Way 
and the host galaxy. 
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the photon timing information within an exposure to reduce telemetry volume. At some 
phases of the observation both modes were operated simultaneously through different sized 
spatial windows that were selected in combination with the spatial binning of the image mode 
data so as to match the overall data volume with the available telemetry rate. Event mode, 
and larger spatial windows, are emphasised at the beginning of the observing sequence, when 
the afterglow brightness is changing rapidly and its position is less well known. Later in the 
sequence only image mode data are taken. The intrinsic spatial pixel size of the detector is 
approximately 0.5 arc seconds on the sky, but some image mode data were taken with the 
data binned 2x2 to give 1 arc second pixels. 

Examination of the UVOT finding chart exposure reveals a bright new fading source 
within the XRT positional error circle of the burst at RA = 18h 32m 32.62s Dec = +26d 
20' 21.6" (J2000) with an estimated uncertainty of 0.2". This close to the ROTSE position 
given by Rykoff et al. (2005). 

In Figure 4 we show the light curves in each of the UVOT bands as a function of time. 
These data were derived from measurements made in a 6 arcsec aperture, with background 
determined from a nearby source-free region. Since the UVOT detector counts individual 
photons, it is subject to coincidence loss (pile- up), which is noticeable for count rates above 
10 per second. We have applied the standard coicidence loss correction to the data, derived 
from ground calibrations (Breeveld et al. 2005). The data have also been corrected for 
Galactic extinction, adopting a value of E(B-V)=0.095 for this field (Schlegel, Finkbeiner 
& Davis 1998). We used the data of Schlegel et al. (1998) and Pei et al. (1992) to derive 
extinction values for each filter. Expressed as magnitudes these are 0.688 (UVW2), 0.926 
(UVM2), 0.740 (UVW1), 0.518 (U), 0.412 (B) & 0.316 (V). The data are tabulated in Table 3, 
separated by filter. The detector count rates have been converted to magnitudes and fluxes 
based on in-orbit calibration data. Systematic uncertainties between the filters are estimated 
to be better than 0.15 magnitudes. 

After the source had faded below the threshold of detectability through the UVOT 
color filters, we made observations in unfiltered (white) light. The white light data provided 
detections of the afterglow beyond 10 5 s. We place the white light data on the V magnitude 
and flux scale using a color dependant correction that is derived by relating the measured 
ratio of V to white counts of field stars and photometric standards to their apparent B-V 
color (i.e. B-V color not corrected for reddening). The white light data are contained in 
Table 4. 
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2.3.1. Temporal Analysis 

We fit a single power law to the light curve for each UVOT broadband filter individually, 
initially using only the data at t < 1000 s after the BAT trigger where the sampling for each 
individual filter is relatively dense. Each filter can be reasonably fit to a common slope of 
a = — 1.14±0.01. We conclude that there is no evidence for a gross change in the optical/UV 
spectrum associated with the temporal decay. 

To maximise the light curve sampling, we combine the data from all the UVOT filters, 
multiplying the individual curves by a factor that normalises them to a common flux scale. 
The results are shown in Figure 5. We also include for comparison the X-ray decay curve of 
Figure 3, which has been multiplied by a constant factor in Figure 5 for display purposes. It 
can be seen that the X-ray and optical/UV signals initially fade at a similar rate, but that 
the optical/UV flux flattens compared to the X-ray curve in the interval 1000-10000s (i.e. 
the ratio of optical to X-ray flux increases with time). After about 10000s the optical curve 
steepens in a similar way to the X-ray curve. 

The complex nature of the optical/UV light curve is confirmed by the fact that a single 
power law does not provide a good fit to the combined UVOT data over the entire time 
range. The chi-square for the fit is 1038 for 115 degrees of freedom. We instead tried various 
combinations of power-laws and found that the best fit was obtained when we fit three 
power law components simultaneously. The model consists of the sum of two power laws 
(Holland et al. 2004) that breaks to a third, single power law at late time. This is physically 
consistent with the flux being dominated by a reverse shock at very early time, a forward 
shock at intermediate time, and a jet break at late time. The best fitting model had an 
early-time power law decay of aii = — 1.56^08 transitioning to a power law with a slope of 
a 2 = —0.62^04 at t t = A32tH 3 s. This transition time is the time at which both power 
law components contribute equally to the observed flux from the afterglow. An achromatic 
break occurs at tj = 13 133^gg} s after the BAT trigger. We interpret this as a jet break. 
The decay index after the jet break is aj = — 1.76iooi- The time of the jet break, and 
the post-break index, are consistent with the values from the X-ray data. The chi-square 
value for this model is 555 for 111 degrees of freedom, the high value reflecting significant 
short-term scatter of individual points about the model. This model is compared to the 
individual filter data in Figure 4. 

Klotz et al. (2005) find that the i?-band light curve of the GRB 050525A afterglow 
underwent a sudden rebrightening of 0.65 mag at 1968 s after the BAT trigger. This was 
at a time when GRB050525a was not visible to Swift. We show the Klotz et al. data in 
Figure 6, normalised to the UVOT curve. If we add this rebrightening to our model the best 
fit has en = -2.14lg;|, a 2 = -1.04+°;^, aj = -1.73^, t t = 39+^ s, and tj = 17747^ 
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s with a chi-square of 597 for 111 degrees of freedom. Although the chi-square value is 
formally larger with a rebrightening, an examination of Figures 4 & 6 suggests that both 
models provide comparable fits. 

Soderberg et al. (2005) have reported a late-time measurement of the GRB050525a 
afterglow brightness using HST/ACS. This is illustrated in Figure 6. The HST measurement 
is approximately 3 magnitudes brighter than predicted at this time by an extrapolation of 
our best fit afterglow model. If the HST measurement represents a constant background flux, 
for example due to the host galaxy or slowly varying supernova emission, then we estimate 
that this would contribute about 25% of the flux that we measure in our UVOT white light 
detection at ~ 2.5 x 10 5 s. The effect of such a constant additional flux on the best fit model 
is illustrated in Figure 5 as the dotted line. 



2.3.2. Multi-wavelength Spectral Analysis 

Combining data across the Swift instruments provides a powerful set of diagnostcs for 
prescribing the instantaneous spectral properties of a source and its temporal evolution. In 
order to fit an instantaneous spectral model to the UVOT data, it is necessary to interpolate 
count rates through each filter to a common epoch. We build a broad-band spectrum of 
the source at T+250s using the appropriate decay index for each filter, and the nearest 
exposure to T+250s as a normalization reference. The XRT spectrum described in section 
2.2.2 was also renormalized to T+250s according to the early decay slope from Figure 3. 
As prior assumptions within the spectral model we take the XRT slope and 2-10 keV flux 
determined in section 2.2.2, the Galactic Hydrogen column density in the direction of the 
source of 9 x 10 20 cm~ 2 (Dickey & Lockman 1990) and the Galactic extinction in the same 
direction characterized by Pei (1992) with E(B-V) = 0.095 (Schlegel et al 1988) and R v = 3.1. 
Applying this model directly to the combined UVOT and XRT spectrum at T+250s yields 
a poor fit with x 2 — 5.3 x 10 4 for 24 dof. This particular model systematically over-predicts 
the optical and UV flux. We attempt to correct this by first adding dust extinction and 
gas absorption from the host galaxy at a fixed redshift of z — 0.606. The dust is initially 
assumed to have a content identical to the SMC with Ry = 2.93 (Pei 1992). This does not 
provide an acceptable fit, with x 2 — 272 for 22 dof, where Nn(host) = (1.8 ± 0.8) x 10 21 
cm~ 2 and E(B-V)(host) = 0.43 ± 0.03. Adopting host dust characteristics identical to the 
Galaxy and LMC populations (Pei 1992) does not improve the fit significantly, yielding x 2 
= 208 and 242, respectively for 22 dof. An acceptable fit is, however, obtained if we add a 
spectral break to the model at 1 keV. Details of the fits to the model with a spectral break 
are provided in Table 5. The table includes results obtained using SMC, LMC and Galaxy 
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dust distribitions. The SMC prescription is preferred statistically over the other two. 

We then perform spectral fits using the broken powerlaw model with SMC-like dust in 
the host galaxy for two other epochs, T+800s and T+25000s. Data and best-fit models are 
plotted in Figure 7 while best-fit parameters are listed in Table 6. Other than temporal 
decay, we find no evidence for spectral evolution between 200s and 800s after the burst; the 
spectral break has not moved within errors. However at T+25,000s the spectrum no longer 
requires a break between the XRT and UVOT energy bands. The UVOT spectral index 
is identical to the XRT index, within measurement uncertainties. There is no evidence for 
evolution in the host's gas or dust enviroments over these epochs. 

3. Discussion 

GRB050525a is the second most fluent GRB to have been observed by Swift to date, and 
is the first bright low redshift burst to have been observed since all three Swift instruments 
have been operational. The hard X-ray /Gamma-ray temporal charateristics and spectrum 
of this bright burst have been well characterised by the BAT. Added to this, the XRT and 
UVOT instruments provide well sampled temporal multiwavelength decay light curves and 
spectra starting about 1 minute after the burst trigger and extending to several days after 
the burst, making this one of the best covered GRBs thus far. The ratio of the gamma-ray 
fluence to the X-ray flux is typical of the general population of long bursts observed in the 
pre-Swift era (cf. Roming et al. 2005b). 

In an initial period of about 4 hours after the burst, the decline of the X-ray afterglow 
flux of GRB050525a can be well represented by a simple power law with an index —1.2, 
apart from a probable flare which occurs after about 300s. The optical temporal signature 
is more complex, however, and exhibits a significant flattening over that same period. It 
can be reproduced by the combination of two power laws (or two power laws combined with 
a re-brightening episode), one declining more rapidly than the X-ray flux, and the other 
more slowly. This can be interpreted as the combintion of a steep reverse shock component, 
combined with the flatter optical apparition of the forward shock seen in X-rays. The X- 
ray and optical/UV multiwavelength spectral fits suggest that the cooling frequency of the 
electrons in the forward shock has moved through the optical band by 25000s after the burst. 
There is good evidence for a break in the time evolution of both the X-ray and optical/UV 
data at 13000-14000s (~0.15 days) after the trigger. The break is consistent with being 
achromatic i.e. the break time is, within errors, consistent between the X-ray and optical 
bands, and the post-break slope is also consistent within the errors (cf. sections 2.2.1 and 
2.3.1). Its achromatic nature suggests that this is the so-called jet break, which occurs when 
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Table 4. UVOT White Light 



T(start) 


T(stop) 


Exposure 


V Mag a 


V Flux b 




(seconds) 








93167 


111189 


2378 


21.35 ± 0.27 


0.103 ± 0.025 


173859 


329532 


11750 


22.82 ± 0.43 


0.027 ± 0.010 


329532 


2599088 


146267 


> 24.11 


< 0.008 



Equivalent V-band magnitude 

Equivalent V-band flux in units of 10~ 16 erg cm -2 s _1 A 



Table 5. SED at T+250s 



SMC LMC Galaxy 



JV H (host) (cm" 2 ) 1.2±0.8xl0 21 1.2±0.8xl0 21 1.2±0.8xl0 2 

Pi -0.60 ± 0.04 -0.63 ± 0.03 -0.62 ± 0.04 

/? 2 * -0.97 -0.97 -0.97 

Ebroak (keV) 1.0 ±0.2 1.0 ±0.2 1.0 + / -0.2 

E(B-V)(host) 0.04 ± 0.02 0.06 ± 0.03 0.06 ± 0.04 

X 2 (dof) 27(20) 31(20) 41(20) 



* Fixed parameter 

Note. - Best fit spectral parameters for the Spectral Energy Distribution 
(SED) of GRB 050525a at T+250s for three dust populations assumed for the 
host galaxy. The model continuum is a broken powerlaw {f3\ where E < Ebroak 
and fi 2 where E > E broak ). The slope in the X-ray band, /3 2 = 0.97, and the 
2-10 keV flux are fixed at the XRT determined values. Galactic values for gas 
absorption and dust extinction are used, while the gas and dust properties of the 
host galaxy (Nn(host) and E(B-V)(host)) are allowed to float. 
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the beaming angle of the decelerating relativistic flow approaches the colimation angle of the 
jet (Sari, Piran & Halpern 1999). 

Comparing the temporal and spectral indices in both the UVOT band and the XRT 
band with the theoretical models (e.g. Table 1 in Zhang & Meszaros 2004) as well as the 
ordering of various temporal segments on the X-ray and optical lightcurves (Sari et al. 1998; 
Chevalier & Li 2000), we find that the global data are not consistent with a model that invokes 
a massive stellar wind. This is also supported by the fact that the spectral break energy 
Eb rea k in our joint XRT-UVOT spectrum (which is consistent with a cooling break) evolves 
downwards with time, since v c increases with time in the wind model. The UVOT-XRT data 
before ~ 400s could be accommodated with a slow-cooling wind model in which the UVOT 
band is between v m and v c (with the temporal index (1 — 3p)/4) while the X-ray band is above 
v c (with the temporal index (2 — 3p)/4). However, the v m < v < u c segment is expected to 
happen at a later time, and it is typically the last segment in the lightcurve so that it is very 
unlikely that there would be a switch back to a shallower decay in the wind model (Fig. 1 of 
Chevalier & Li 2000). The early steep-to-shallow transition observed in the UVOT lightcurve 
naturally rules out the wind interpretation, and is consistent with a reverse-forward shock 
model as discussed below. Another possibility is that the fireball may be initially in a wind 
medium, but later runs into an ISM after it crosses a wind-termination shock. An additional 
motivation for this interpretation is the optical-rebrightening at ~ 2000s revealed by the 
TAROT data (Klotz et al. 2005), which could be presumably the signature of the termination 
shock. Although the possibility that the rebrightening is due to a density enhancement in 
the medium surrounding the burst is promising, the above scenario lacks direct proof since 
the wind signature, if it is one, terminates at a much earlier time (e.g. ~ 400s) than the 
epoch of the optical rebrightening bump. 

Instead, the XRT and UVOT data before the presumed jet break at 0.15 d are in good 
agreement with the standard fireball model for a constant ISM density (e.g. Sari, Piran 
& Narayan 1998). According to this model, the spectral indices {(3) above and below the 
synchrotron cooling frequency v c are — p/2 and — (p — l)/2, respectively; while the temporal 
indices (a) for a band which is above and below v c are (2 — 3p)/4 and (3 — 3p)/4, respectively. 
Interpreting E break as the cooling break, and taking p = 2.2, one expects ax = —1.15, 
a = —0.9, p x = —1-1, and (3 — —0.6. All these are in general agreement with the data, 
particularly when the optical re-brightening at ~ 2000s (Klotz et al. 2005) is taken into 
consideration when fitting model parameters. At T + 25, 000s, the energy spectrum of the 
burst is consistent with a single power law spectrum extending between the XRT and UVOT 
bands, with index f3 ~ —0.95. This is consistent with the fact that the cooling frequency has 
already crossed the UVOT band at this epoch. The decay slope expected in the UVOT band 
at this time should be ~ —1.15 in the standard model. The observed decay slope, however, 
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is steeper (~ —1.6). This is due to the jet break around ~ T + 14000s. According to the 
standard afterglow model, v c evolves with time as oc t~ x l 2 . The data, on the other hand, 
requires a slightly faster evolution of E break . Enhanced cooling is thus needed, possibly with 
an evolving e#, the magnetic equipartition parameter in the shock. An evolving e# was also 
invoked to interpret data GRB 050128 (Campana et al. 2005). 

A temporal break is identified in the X-ray lightcurve at around 14000s (0.15d). The 
UVOT data show a break around the same time. Since there is no apparent spectral change 
aross the break, and since the break time is "achromatic", this points to a jet break, corre- 
sponding to the time when the fireball Lorentz factor 7 ~ l/0j (Rhoads 1999; Sari, Piran 
& Halpern 1999). The post-break temporal index, however, is shallow (~ —1.6) compared 
with the model prediction (which is ~ — p = —2.2). Such a shallow jet break has also been 
seen previously in some other bursts (e.g. GRB 990123, Kulkarni et al. 1999). It might 
be shallow because the sideways expansion effect is not significant (Panaitescu & Meszaros 
1999). Another possibility is that the jet break turn-over time may be finite (Panaitescu & 
Meszaros 1999), so that the post-break asymptotic regime has not yet been reached. If this 
is the case, then the effective jet break time could be later than we have fitted. For example, 
while the increased number of free parameters are not justified by the quality of our data, we 
could force a 3-power-law fit to our X-ray data with a second break to a slope of a — —2.2 
at about T= 1 — 1.5 x 10 5 s, yielding an 'effective' single break time of about 50, 000-60, 000s. 
This is illustrated in Figure 3, where the dash-dot line shows a break to a slope of a — —2.2 
at T= 1.5 x 10 5 s. The final two XRT points in isolation are equally consistent with a slope 
of a = —1.6 and a = —2.2. It is also still possible that the observed temporal break is due 
to reasons other than a jet break (e.g. Gendre & Boer 2005), but given that this is the first 
time one sees an achromatic break in both the X-ray and the optical band, we tentatively 
conclude that the current data support a jet break interpretation. 

If the break is indeed attributed to a jet, one can derive a jet angle (e.g., Sari et al. 
1999) 

B 3 = 0.12 ( JL.) ^ ( Wl.) V8 = 5.6 x IP" 2 (^) ^ ( \ ) "* , (1) 
3 \1 + z) \E iso , 53 J V 0.6 7 V0.15da y ; ' yj 

where tj is expressed in days, no is the ambient density, w 7 is the ejecta-to-gamma-ray 
efficiency, and £"180,53 is the isotropic energy in units of 10 53 erg. In our calculations we use 
tiq = 3 cm -3 and r] 1 = 0.2 (Ghirlanda, Ghisellini & Lazzati 2004). Therefore 6j ~ 3.2° 
if we use a jet break time of tj ~ 0.15 day, or 9j ~ 5.4° if we adopt a more conservative 
jet break time of tj ~ 0.6 day corresponding to a finite roll-over as discussed above; these 
two values bracket the possible range of tj. This angle could be interpreted as the physical 
opening angle of a jet if the jet has a uniform distribution of energy, or as the observer's 
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viewing angle with respect to the jet axis in a structured jet model (Rossi et al. 2002; Zhang 
& Meszaros 2002; Kumar & Granot 2003). Within the uniform jet model, we estimate 
the actual gamma-ray energy emitted to be E 1 = E ison (l — cos9j) = 3.6 x 10 49 ergs for 
tj rsj 0.15 day or = 1.0 x 10 50 ergs for tj ~ 0.6 day. 

We can combine the information on the jet with the measurements of the redshift, 
E p , E is0tl , tj and the lag between the high energy and low energy gamma-ray light curves 
(section 2.1) to test whether this burst is consistent with various proposed empirical re- 
lations. Four relations involve E pt , the peak energy in the burst frame for the entire 
burst: the Amati relation — -Ept.Amati oc E^ (Amati et al. 2002); the Ghirlanda relation — 
£ P t,Ghirianda oc E^ 7 (Ghirlanda et al. 2004); the Yonetoku relation — -E pt ,Yonctoku oc L9£ 
(Yonetoku et al. 2004); and the Liang-Zhang relation — -E pt ,Lz oc E^ so t Q j M . (Liang & Zhang 
2005). For GRB 050525a's values of E iso , E 7 , tj and L iso , we calculate £pt,Amati=144 keV, 
£ P t,Ghirianda=46 keV (if tj = 0.15 day 96 keV if tj = 0.6 day), -E pt ,Yonctoku = 186 keV, and 
-£'pt,LZ=50 keV (if tj = 0.15 day, 122 keV if tj = 0.6 day), which should be compared to 
the observed value of E pt = 126.6 ± 5.5 keV. The set of bursts used to calibrate these rela- 
tions is small enough that each additional burst with a redshift and a spectrum will change 
the relations quantitatively; in addition, there is significant scatter around these relations 
(G. Ghirlanda, 2005, personal communication). 

According to the lag-luminosity relation, shorter lags between the emission in two 
gamma-ray bands are correlated with larger L iso (Norris, Marani & Bonnell 2000). Using the 
methodology in Band, Norris & Bonnell (2004), we calculate a redshift z = 0.69 ± 0.02 for 
GRB050525a from the observed values of the peak flux, lag and spectrum at the peak of the 
lightcurve (Norris et al. 2005); the uncertainty in this derived time-lag redshift accounts for 
the uncertainty in the lag and the peak flux, but not for the systematic uncertainty resulting 
from the somewhat different energy bands used to calibrate the lag-luminosity relation and 
for the Swift observation. Given the large redshift range over which bursts are detected, the 
agreement between this time-lag redshift and the spectroscopic redshift is impressive. 

The UVOT lightcurves are consistent with the existence of an early reverse shock compo- 
nent (Meszaros & Rees 1997; Sari & Piran 1999). The best-fit initial decay index is sensitive 
to whether we include a re-brightening episode in our model, as suggested by Klotz et al. 
(2005). It ranges from c*i ~ —1.5 without re-brightening to c*i ~ —2.1 with re-brightening. 
This is in the range observed for previous suggested detections of the reverse shock in GRB 
990123 (Akerlof 1999) and GRB 021211 (Fox et al. 2003; Li et al. 2003), which are typically 
~ —1.9. A shallower decay corresponds to a shallower electron index in the reverse shock 
region, p = (— 4c* i — l)/3 ~ 1.75 (Zhang, Kobayashi & Meszaros 2003). In contrast with 
Shao & Dai (2005) who interpreted the tentative jet break as the forward shock peak, in 
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our best fit, the forward shock peaks at a much earlier time, before 400 s. This corresponds 
to a typical "flattening" -type early afterglow (Zhang et al. 2003), which usually requires a 
magnetized central engine (see also Fan et al. 2002; Kumar & Panaitescu 2003). The best fit 
temporal index that we derive for the forward shock component in the optical/UV without 
a re-brightening is -0.62, flatter than the expectation of the simplest model (-0.9). However 
if we include a re-brightening episode as described above, the slope of the forward shock 
component steepens to —1.04, which would be more consistent with the simple model. 

The Swift programme is supported by NASA, PPARC and ASI (contract number 
I/R/039/04). 
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Table 6. The SMC-dust model from Table 5 applied to the GRB 050525a SED at three 

epochs. 



T+250s T+800s T+25,000s 



N H (host) (cm- 2 ) 1.2±0.8xl0 21 1.5±0.9xl0 21 1.8±1.0xl0 2 

(3i -0.60 ±0.04 -0.67 ± 0.06 -0.94 ±0.10 

(3 2 * -0.97 -0.97 -0.97 

Ebroak (keV) 1.0 ±0.2 0.9 ± 0.2 

E(B-V)(host) 0.04 ± 0.02 0.09 ± 0.04 0.08 ± 0.05 

X 2 (dof) 27(20) 18(20) 15(20) 



* Fixed parameter 
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